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We report the first measurement of the net-charge fluctuations in Pb-Pb collisions at -^snn = 
2.76 TeV, measured with the ALICE detector at the CERN Large Hadron Collider. The dynamical 
fluctuations per unit entropy are observed to decrease when going from peripheral to central colli- 
sions. An additional reduction in the amount of fluctuations is seen in comparison to the results 
from lower energies. We examine the dependence of fluctuations on the pseudo-rapidity interval, 
which may account for the dilution of fluctuations during the evolution of the system. We find 
that the ALICE data points are between the theoretically predicted values for a hadron gas and a 
Quark-Gluon Plasma. 



PACS numbers: 25.75.-q,25.75.Nq,12.38.Mh 



The ALICE experiment [lj at the Large Hadron Col- lider (LHC) is a multi-purpose detector designed to study 
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the formation and evolution of nuclear matter at high 
temperatures and energy densities. One of the major 
goals of the experiment is to explore as many signals 
as possible towards characterizing the properties of the 
Quark-Gluon Plasma (QGP), the deconfined state of 
quarks and gluons, produced in high energy heavy-ion 
collisions. The study of event-by-event fluctuations pro- 
vides a powerful tool to characterize the thermodynamic 
properties of the system. The fluctuations of conserved 
quantities, like net-charge of the system, are predicted to 
be one of the most sensitive signals of the QGP formation 
and phase transition, and may provide complementary 
understanding of strong interactions ■ 

In the QGP phase, the charge carriers are quarks with 
fractional charges, whereas the particles in a hadron gas 
(HG) carry unit charge. The fluctuations in the net- 
charge depend on the squares of the charge states present 
in the system. Consequently, the net-charge fluctuations 
in the QGP phase are significantly smaller compared to 
that of a HG Q. At the same time, if the initial QGP 
phase is strongly gluon dominated, the fluctuation per 
entropy may further be reduced as the hadronization of 
gluons increases the entropy Thus the net-charge 
fluctuations are strongly dependent on which phase they 
originate from. However, the net-charge fluctuations 
may get affected by uncertainties arising from volume 
fluctuations, so one considers the fluctuations of the ratio, 
R = N + /N^. Here N + and iV_ are the numbers of pos- 
itive and negative particles respectively, measured in a 
specific transverse momentum (pr) an d pseudo-rapidity 
(77) window. The parameter R is related to the fluctu- 
ations of the net-charge via the £>-measure as per the 
following expression 0, 0, [ft : 



D= (N ch )(6R 2 



(SQ 2 



(1) 



which provides a measure of the charge fluctuations per 
unit entropy. Here the (...) denotes an average of the 
quantity over an ensemble of events. The term (8Q 2 ) is 
the variance of net charge, Q = N + — N- and iV c h = 
N+ + N_ . The D-measure has been estimated for sev- 
eral different theoretical considerations including those of 
the lattice calculations. In a simple picture by neglect- 
ing quark-quark interactions, D is found to be approxi- 
mately 4 times smaller for a QGP compared to a HG Q. 
Lattice calculations which include the quark-quark in- 
teractions give a quantitatively different estimate for a 
QGP phase, still significantly smaller than for a HG. It 
has been shown that D = 4 for an uncorrelated pion gas, 
and after taking resonance yields into account, the value 
decreases to D ~ 3. For a QGP, D is significantly lower 
and has been calculated to be D ~ 1.0-1.5 where the 
uncertainty arises from the uncertainty of relating the 
entropy to the number of charged particles in the final 
state [5J. Thus, a measurement of D can be effectively 
used as a probe for distinguishing the two phases, the HG 



and the QGP. However in reality, these fluctuations may 
get diluted in the rapidly expanding medium due to dif- 
fusion of particles in rapidity space @, [t| . Several other 
effects, such as collision dynamics, radial flow and final 
state interactions may also affect the amount of measured 
fluctuations 0, [l(J 11 1 . 



In the exp eriment, the net-charge fluctuations are best 



studied [llHK 
fined as: 

f(+-,dyn) 



by calculating the quantity Vr_\ dyn) de 

(N + (N+-1)) , (iV_(iV_ - 1)) 



(iV 4 



(AM 2 
(N-N + ) 



"{N-){N+y (2) 

which is a measure of the relative correlation strength 
of particle pairs. A negative value of vt_\ ^yn) signi- 
fies the dominant contribution from correlations between 
pairs of opposite charges. On the other hand, a positive 
value indicates the significance of the same charge pair 

correlations. The Vr^ i( jyn) has been found to be robust 

against random efficiency losses [16l4l8| . Z?-measure and 
dyn ) are related to each other by [5[: 

(A r d>(+-,dyn) « D - 4. (3) 

The values of djfa \ need to be corrected for global 

charge conservation [la]. The predictions for the D- 
measure are based on the assumption of vanishing net- 
charge in the system. However, in a realistic situation, 
the system under consideration has a small but finite net- 
charge. A correction due to finite net-charge effect also 
needs to be applied 0]. 

In this letter, we report the first measurements of net- 
charge fluctuations, by calculating vr_\ dyn) and the D- 

measure, as a function of collision centrality in Pb-Pb 
collisions at ^/snn =2.76 TcV at the LHC with the AL- 
ICE detector. We also make a comparison of the experi- 
mental results to the theoretical predictions. 

Details of the ALICE experiment and its detectors can 
be found in [l|. For this analysis, we have used the Time 
Projection Chamber (TPC) [lj| to reconstruct charged 
particle tracks. The detector provides a uniform accep- 
tance with an almost constant tracking efficiency of about 
80% in the analyzed phase space (\r]\ < 0.8 and 0.2 GeV/c 
< Pt < 5 GcV/c). The interaction vertex was measured 
using the Silicon Pixel Detector (SPD), the innermost 
detector of the Inner Tracking System (ITS) of ALICE. 
In the analysis, we have considered events with a vertex 
\v z \ < 10 cm to ensure a uniform acceptance in the cen- 
tral pseudo-rapidity region. The minimum bias trigger 
of ALICE consisted of a coincidence of at least one hit on 
each of the two VZERO scintillator detectors, positioned 
on both sides of the interaction point, while at the startup 
of data taking period an additional requirement of hav- 
ing a coincidence with a signal from the SPD was also in- 
troduced. The background events coming from parasitic 
beam interactions arc removed by a standard offline event 
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FIG. 1. Dynamical net-charge fluctuations, d yn ) and 

their corrected values, dyn)j f° r charged particles pro- 

duced in Pb-Pb collisions at ^/snn = 2.76 TeV as a function 
of centrality expressed in terms of the number of participating 
nucleons. rfTL dyn ) points are shifted along x-axis for better 
representation. Superimposed are the results for pp collisions 
at y/s — 2.76 TeV. The results are shown for two different 
Arj windows. The statistical (bar) and systematic (box) er- 
rors are plotted. 



calculate the i( jyn) from the experimental measure- 
ments of positive and negative charged particles counted 
in A?y windows, defined around mid-rapidity (for ex- 
ample, A?y = 1 corresponds to —0.5 < r\ < 0.5) and 
in the px range from 0.2 to 5.0 GeV/c. In Figure [TJ 
we present the dyn) as a function of centrality, ex- 
pressed in terms of the number of participating nucleons. 
Moving from left to right along the x-axis of the figure 
corresponds to moving from peripheral to central colli- 
sions. The results are presented for Arj = 1 and 1.6, for 
both Pb-Pb and pp collisions. In all cases, the magni- 
tude of z^(_j dyn) is observed to be negative, indicating 

the dominance of the correlation term in Eq. 2. The 

absolute values of d yn) f° r PP collisions are larger 

compared to those measured for Pb-Pb collisions. When 
going from peripheral to central events, the absolute val- 
ues of dyn) are seen to decrease monotonically. 

The values of ^( + _ idyn ) have to be corrected for global 
charge conservation and finite acceptance [lij]. If all 
charges are accepted, the global charge conservation 
would lead to vanishing fluctuations. This will yield 

the minimum value of dyn) to be -4/(A r to tai), where 

(^Vtotai) is the average total number of charged particles 
multiplicity produced over full phase space. The cor- 



rected value of v. 



(+- ,dyn) 



is then: 



selection procedure, which requires the VZERO timing 
information and hits in the SPD. 

We present the results as a function of centrality that 
reflects the collision geometry. The collision centrality 
is determined by cuts on the VZERO multiplicity as de- 
scribed in [2^ . A study based on Glauber model fits plT - 
231 ] to the multiplicity distribution in the region corre- 
sponding to 90% of the most central collisions, where the 
vertex reconstruction is fully efficient, facilitates the de- 
termination of the cross section percentile and the num- 
ber of participants. The resolution in centrality is found 
to be < 0.5% RMS for the most central (0-5%) collisions, 
increasing towards 2% RMS for peripheral (70-80%) col- 
lisions [2Cf . 

We require tracks in the TPC to have at least 80 re- 
constructed space points with a \ 2 P er TPC cluster of 
the momentum fit less than 4. We reject tracks with dis- 
tance of closest approach (dca) to the vertex larger than 
3 cm both in the transverse plane and in the longitudi- 
nal direction. We have performed an alternative analysis 
with tracks reconstructed using the combined tracking of 
ITS and TPC. In this case, the dca cuts were 0.3 cm in 
the transverse plane as well as in the longitudinal direc- 
tion. The results obtained with both tracking approaches 
are in agreement. 

The data analysis has been performed for Pb-Pb colli- 
sions at y / 5nn = 2.76 TeV and pp collisions at the same 
centre-of-mass energy. An identical analysis procedure 
has been followed for both Pb-Pb and pp data. We 



,,corr 

v {-{ — ,dyn) — v {+- ,dyn) 



total/ 



(4) 



Since the value of (-ZVtotai)j has not been reported by 
experiments, we have obtained these values from HI- 
JING Q and PYTHIA [H] event generators for Pb-Pb 
and pp collisions, respectively. The corrected values, 



'(+- ,dyn): 



are calculated and plotted in Figure Q] as a 
function of the number of participating nucleons for 
Pb-Pb as well as pp collisions. The absolute values of 

dyn) are smaller compared to .dyn) in all cases. 

The differences are more apparent for pp and peripheral 
Pb-Pb collisions than for central collisions. 

The calculations of the D-measure were consequently 

done starting from the corrected values of dyn ) in 

Eq. 3. A systematic check of correcting the D-measure 
has been performed by using D— 7>_D/(C M C I? ), where — 



(N. 



±4- and C r , 



1 



- $g^y 0- A difference of 3-7% (de- 
pending on the Ary window) has been included as one of 
the systematic errors to D. For the rest of the manuscript 

the corrected values of dyn ) and D are presented. 

The systematic uncertainties have additional contri- 
butions from the following sources: (a) uncertainty in 
the determination of the interaction vertex, (b) differ- 
ent magnetic field polarities, (c) contamination from sec- 
ondary tracks {dca cuts), (d) centrality definition using 
different detectors, (e) selection criteria at the track level, 
and (f) different tracking scenarios. The total systematic 



error on v 



amounts to 6-10% in going from pe- 
ripheral to central collisions. The error on the product of 



(+- ,dyn) 
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FIG. 2. 



(+-,dyn) 



(left axis) and D (right axis) as a 



function of the number of participants for Arj = 1 and 1.6 
in Pb-Pb collisions at ^/snn = 2.76 TeV and pp collisions 
at y/a — 2.76 TeV. Also shown are results from the HIJING 
event generator for both the A77 windows. Both statistical 
(error bar) and systematic (box) errors are shown. 
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FIG. 3. {Nch}v?+l idyil ) (left axis) and D (right axis) as a 
function of Ar] window for three different centrality bins in 
Pb-Pb collisions at ^/snn = 2.76 TeV. The data points are fit- 
ted with the functional form, erf(A?j/v8o'f)- The dashed lines 
correspond to the extrapolation of the fitted curves. Both sta- 
tistical (error bar) and systematic (box) errors are shown. 



number of charged particles and 1^°" dyn ^ remains within 
7-13% at all centralitics. The systematic and statistical 
uncertainties in all the figures are represented by boxes 
and error bars, respectively. The statistical errors are 
small and within the sizes of the symbols in most cases. 



Figure [5] presents the values of (iVch) 



'(+-,dyn) 



and D 



in the left and right axes, respectively, as a function of 
the number of participating nucleons. The (iV c h) values 
have been measured for different ccntralities and A77 win- 
dows, and corrected for detector inefficiencies [2(|. Both 
the results from the Pb-Pb and pp analyses are shown. 
The shaded bands in the figure indicate the predictions 
for a HG and a QGP. The figure also shows the results 
from the HIJING event generator for A77 = 1 and 1.6, 
which are observed to be close to the HG line and at the 
same time independent of centrality. The pp results for 
Arj = 1.6 agrees very well with the HG prediction. The 
experimental results for Pb-Pb for both the A77 windows 
are observed to be below the HG predictions and above 
those of the QGP. The values of D for A?/ =1.6 are lower 
compared to those for Arj = 1 for all centralitics. 

A decreasing trend of D has been observed while going 
from peripheral to central collisions, as seen in Figure [5] 
This centrality dependence may arise partly because of 
the presence of radial flow [l(| ■ The radial flow velocity 
could lead to the kinetic focussing of the produced par- 
ticles, causing a narrowing of the opening angles. There- 
fore, it is expected that the number of positive and neg- 
ative particles may get redistributed in a finite phase 
space. This may affect the magnitude of net-charge 

fluctuations. The effect of radial flow on d yn ) has 

been estimated by using two different methods. First, 



we have used an afterburner [26j on the HIJING events 
where the particles get a boost in the transverse mo- 
menta because of the radial flow velocity. The magni- 



tudes of v. 



(+-,dyn) 



and D are observed to be close to 



each other for both HIJING and HIJING with the af- 
terburner. In the second method, D is calculated using 
the AMPT model 27J. Both versions of this generator, 
the string melting and the default, were studied. The 
observed centrality dependence in the data are not seen 
with the AMPT model. These studies indicate that the 
presence of radial flow may not be responsible for the 
centrality dependence of the £)-measurc. 

The measured fluctuations may get diluted during the 
evolution of the system from hadronization to kinetic 
freeze-out because of the diffusion of charged hadrons 
in rapidity. In ref. 0, this has been addressed 
where a diffusion equation has been proposed to study 
the rapidity dependence of net-charge fluctuations. It 
has been conjectured that, taking the dissipation into 
account, the asymptotic value of fluctuations may be 
close to the primordial fluctuations in large rapidity win- 
dows corresponding to the hadronization stage. This 
has been explored for the ALICE data points by plot- 



ting (N, 



(+-,dyn) 



and D as a function of A77 for three 
centrality bins, as shown in Figure [3] We observe that 
for a given centrality bin, the .D-measure shows a strong 
decreasing trend with the increase of A77. In fact, the cur- 
vature of D has a decreasing slope with a flattening ten- 
dency at large A77 windows. Following the prescriptions 
of [H ||| , we fit the data points with the functional form, 
crf(Ar]/\/8<jf ), which represents the diffusion in rapidity 
space. Here, at characterizes the diffusion at freeze-out. 
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FIG. 4. Energy dependence of the net-charge fluctuations, 



measured in terms of (iV c h}i 



( + - ,dyn) 



(left axis) and D (right 



axis), for the top central collisions. The results from the 
STAR [llj and ALICE experiments are presented for Ar/ — 1 
after the correction for charge conservation. The ALICE re- 
sults are further extrapolated to obtain the asymptotic value, 
shown by the solid red circle. Both statistical (error bar) and 
systematic (box) errors are shown. 



The resulting values of <r f are 0.42 ± 0.04, 0.48 ± 0.06 and 
0.53 ±0.05 for the 0-5%, 30-40% and 50-60% centralities, 
respectively. The fitted curves are shown as solid lines 
in Figure [3l The dashed lines are extrapolations of the 
fitted curves to higher A77, which yield the asymptotic 
values of D. For the top 5% centrality, the measured 
values of D are 2.47± 0.01(stat.) ± 0.13(sys.) for At? = 1 
and 2.10 ± 0.02(stat.) ± 0.12(sys.) for Arj = 1.6. The 
extrapolated value of D is 1.99 ± 0.09(stat.) ± 0.17(sys.). 

The evolution of the net-charge fluctuations with beam 
energy can be studied by combining the ALICE data for 
Pb-Pb collisions with those of the STAR experiment [ll[ 
for Au-Au collisions at four RHIC energies. In Figure |H 
we present the values of (-/V c h)i^ c °" dyn ^ (left axis) and 
D (right axis) for the top central collisions from AL- 
ICE at v /iii' =2.76 TcV and for STAR, Au-Au col- 
lisions at four different energies. The ALICE data points 
correspond to Arj = 1 and 1.6, whereas for STAR the 
values for Arj = 1 are shown. For the STAR data, 
(dN ch /d V )v c { ™ , dyn) is plotted instead of (N ch )u^Z Ayn) , 
where dN c h/ ' dr\ is approximately equal to (N c h) for A77 = 
1 at central rapidity. The theoretical predictions for a 
hadron gas and a QGP are also indicated in the figure. 

In Figure 01 we observe a monotonic decrease in the 
magnitude of the net-charge fluctuations with increas- 
ing beam energy. For the top RHIC energy of ^/snn = 
200 GeV, the measured data of STAR is close to the pre- 
diction for a hadron gas. Below this energy, the data 
points are above the hadron gas value. We note that the 
STAR results are obtained for Ai] = 1 where the fluctu- 



ations may still be significantly affected by diffusion. At 
2.76 TcV, we observe that the values of D for both the 
At] windows are within the predictions for a hadron gas 
and a QGP. At A?y = 1.6 the value of D is closer to that 
of the QGP prediction. 

In summary, we have presented the first measurements 
of dynamic net-charge fluctuations at the LHC in Pb-Pb 

collisions at ^/snn = 2.76 TeV in terms of !( iyn)> 

and their corrected values, ^ c °" dyn ^ (corrected for charge 
conservation and finite acceptance effect). The results for 
pp collisions at the same center-of-mass energy are found 
to be in agreement with hadron gas prediction. The val- 



ues Of Z/( + _. dyn ) 



and 



( + - dyn) 



are seen to be negative in 



all cases, indicating the dominance of the correlation of 
positive and negative charges. A decrease in fluctuations 
is observed while going from peripheral to central colli- 
sions. The £>-measure, which gives the charge fluctua- 
tions per entropy, is calculated from v?^_ dyn \ and from 
the measured average charged particle multiplicity. A 
decreasing trend of D is observed in going from periph- 
eral to central collisions. Model studies indicate that the 
presence of radial flow may not be the cause of this de- 
crease. The dissipation of signal during the evolution of 
the fireball from the hadronization to freeze-out has been 
estimated by fitting D as a function of the Arj window. 
The extrapolation of the fit function yields the asymp- 
totic value of .D, which is not very different from the mea- 
surement at A77 = 1.6. The beam energy dependence of 
charge fluctuations has been studied by comparing the 
ALICE data with those from the STAR experiment at 
RHIC for Au-Au collisions at four energies. A mono- 
tonic decrease in the value of D, measured at Arj = 1, 
has been observed. The STAR data points at RHIC top 
energy are close to the prediction for a hadron gas. This 
may be due to the fact that the fluctuation may be not 
strong enough to be measured or because of the dilution 
of fluctuation during the evolution process. The ALICE 
data points are below the prediction for a hadron gas and 
above that of the QGP. Moreover, these data points show 
an additional decrease of D at Arj = 1.6. For the top cen- 
tral collisions, the measured value of D turns out to be 
2.10 ± 0.02(stat.) ± 0.12(sys.) which is to be compared 
with theoretical predictions. 
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